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ABSTRACT 


The  puleed  propulsive  performance  of  low -thrust  „ taction  jets,  typical 
of  those  used  for  small  spacecraft  attitude  control,  is  analyzed  and  com¬ 
pared  with  the  results  of  laboratory  experiments.  Five  gases,  hydrogen, 
nitrogen,  ammonia,  Freon  =  i2,  and  Freon- 14,  are  investigated  using  a 
48  to  1  expansion  ratio  nozzle.  The  transient  processes  which  dominate  the 
short-pulse  or  limit-cycle  mode  of  thruster  operation  are  formulated.  These 
relationships  show  good  correlation  with  the  data.  The  apparatus,  proce¬ 
dures,  and  tecliniques  required  to  obtain  accurate  test  results  for  a  low- 
thrust,  dynamic  mode  of  operation  are  described.  Impulse  bit  size,  gas 
consumption,  and  specific  impulse  are  characterized  in  terms  of  thruster 
geometry,  gas  properties,  and  command  pulse  width  to  provide  a  basis  for 
optimum  system  design.  A  simplified  method  for  calculating  dynamic  impulse 
bit  size,  dynamic  gas  consumption,  and  pulsed  specific  impulse  as  a  function 
of  command  pulse  width  ia  developed.  Finally,  the  effective  performance  of 
the  gases  tested  is  evaluated  by  a  technique  which  includes  the  influence  of 
tank  and  propellant  weights,  as  well  as  specific  impulse. 
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NOMENC  L.ATURE 


Aoa*{2/v)i/?72Vc 

(Ata*/2Yc)[2/(Y  +  i)]<V»-i)/2(Y-l) 


acoustic  velocity 


equivalent  orifice -to-nozzle  area  ratio 
A,  / ■>  /2  / . ,j_ i  \2(y-l)  a  j 


nozzle  discharge  coefficient 


thrust  coefficient 


thrust  correlation  factor 


exhaust  velocity  coefficient 


diameter 


thrust  level 


proportionality  constant  in  Newton's  second  law 
specific  impulse 

effective  system  specific  impulse 


total  impulse 


polytropic  exponent 
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NOMENCLATURE  (continued) 


empirical  decay  time  integration  factor 

choked  nozzle  flow  factor 

empirical  rise  time  integration  factor 

Mach  number 
molecular  weight 
Knud  sen  number 

Reynolds  number 

pressure 

universal  gas  constant 

temperature 

volume 

specific  volume 
weight  flow  rate 
weight 

pressure  ratio  transformation  =  j[i  -  (Pc/P^) 

[i  -  (p./p,)]1'2 

ratio  of  specific  heats 
nozzle  expansion  ratio 
efficiency 


NOMENCLATURE  (continued) 


8  =  time 

p  =  density 

v  =  standard  deviation  =  0.  7631  X  10'5  lb-sec  or  working  stress 

t  =  dimensionless  time  factor  =  a£0 


SUBSCRIPTS 


a 

= 

ambient 

c 

= 

chamber  or  command 

d 

s' 

decay 

e 

= 

nozzle  exit 

i 

= 

final 

I 

s 

thermal 

i 

s 

initial 

n 

= 

nozzle 

o 

= 

orifice 

P 

= 

polytropic 

r 

= 

reservoir  or  rise 

s 

= 

steady  state 

t 

ST 

throat 

I.  INTRODUCTION 


The  propulsive  performance  of  variouB  gases  with  a  low  expansion 
ratio  (1.  68  to  1)  nozzle  was  reported  in  Refs.  1-3.  In  this  document, 
propulsive  performance  and  design  relationships  for  cold  gas  jet  reaction 
systems  which  might  be  used  in  small  space  vehicle  attitude -control  applica¬ 
tion  are  analytically  derived  and  are  correlated  with  the  results  of  vacuum 
experiments  using  hydrogen,  nitrogen,  ammonia,  Freon-12,  and  Freon- 14 
in  a  high  expansion  ratio  (48  to  1)  noxzle.  These  relationships,  previously 
unavailable  in  systematic  form,  provide  the  analytical  tools  needed  to  design 
low  thrust,  pulsed,  cold  gas  propulsion  systems,  A  simplified  means  :Qf 
calculating  dynamic  impulse  bit  sine,  gas  consumption,  and  specific  Impulse 
is  developed  for  use  in  preliminary  design. 

A  schematic  diagram  of  the  experimental  apparatus  is  shown  in 
Fig.  i.  Pressure  and  thrust  histories  and  the  -weight  of  consumed  gas  were 
measured  so  that  specific  impulse  could  be  determined  as  a  function  of 
command  pulse  width.  A  complete  dei  cription  of  the  experimental  program 
is  presented  following  a  discussion  of  the  analytical  approach. 
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THRUST  STAND  (DETAIL  1} 
SUPPLY  TANK 
CHECK  VALVE 
PRESSURE  GAGE 
PRESSURE  REGULATOR 
THERMOCOUPLE 
SOLENOID  VALVE 
PRESSURE  TRANSDUCER 
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DUCER  (THRUST  PICKUP) 
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CANTILEVER  BEAM  AND 
PNEUMATIC  FEED 


Fig.  i.  Schematic  of  Experimental  Apparatus 
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II.  ANALYTICAL  APPROACH 


A.  PRESSURE  TRANSIENTS 

The  reaction  jet  performance  is  based  on  a  dynamical  analysis  of  the 
thrust  chamber  pressure  history.  In  general,  the  chamber  pressure  history 
is  composed  of  the  steady-state  and  transient  (startup  and  shutdown)  phases 

illustrated  in  Fig.  2.  In  the  following  analysis,  it  is  shown  that  the  transient 

$ 

phase  causes  performance  nonlinearities.  Design  relationships  for  predicting 

and  optimizing  these  nonlinearities  are  developed.  "  >V 

During  startup,  the  flow  rate  through  the  solenoid  supply  valvb  brifiee 

W  exceeds  that  through  the  discharge  nozzle  W  #  causing  thech  amber 
°  “ 

pressure  P  to  build  up.  The  exact  expression  for  the  compressible  isotropic 
pressure  buildup  was  given  in  Ref.  2  as  .  ?  -  ^ 


For  steady-state  operation  dP_/d9  =  0;  and  P  /Py^'P  /Pf  -This  ratio  is  a 
design  variable  (typically  on  the  order  of  0.  99)  and  is  determined  by  the 
selection  of  Aq  and  A{.  '  - 

Equatic  (i)  cannot  be  solved  direedy  but  may  be  simplified  using 

*  i.. 

one  of  three  assumptions:  (i)  the  thruster  nozzle  flow  is  zero  during  startup; 

-v  _ 

(2)  both  nozzles  are  always  choked;  (3)  the  flow  through  the  orifice  oq^Urs 

-  *  C  .  *  ;  - 

at  constant  density.  For  the  case  where  the  thruster  nozzle  flow  is  assumed 


5  *• 


<Vi*  W-WV-K 
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to  be  negligible  during  startup  (assumption  i),  Eq.  (1)  becomes 


ml/py/v  /p„\<Y+1,/Y 


Integrating  Eq.  (2)  between  .the  appropriate  limits  results  in 


P  /P  =  1 
c  s 


[.■IH-rTf 

•» 


2-iv/(v-i) 


(V  -  l)/\  j  '  Bt 


If  both  nozzles  are  assumed  always  choked  (assumption  2),  then  Eq.  (i) 


becomes 


dPc/ae  =  2.2(AoPr/At .  pc)  , 


•  (4) 


which  integrates  to 


Vp»  =  VA« +  [‘  -  y\  -  <VAt>k2T 


Assuming  constant  density  orifice  flow  (assumption  3),  Eq.  (i)  becomes 


dPc/de  =  2.,Rr[i  ■  (Pc/Pr)]: 


r  2  a,P 
2  c 


5 


Equation  (6)  can  be  integrated  directly  to  obtain 


-2t  =  ln{-  y£’  -  By  +  i)  + 


_ 2B _ -2y  -  B  -  (B2  -t  4) 1/2  |  y 

(B2  +  4)1/2  (-y2  -  By  +  1)1/2  I  yr 


To  solve  Eq.  (7)  explicitly  for  P  /F  ,  it  is  assumed  ^that  B^  is  >4,  so 

c  s 

Fc/Ps  =  1  -  (yo  +  B)2  e”2x  +  2B(yo  +  B)e“T  -  B2  (8) 

A  comparison  of  the  results  obtained  with  Eqs.  (3),  (5),  and  (8)  and  meas¬ 
ured  data  indicated  better  correlation  over  a  wider  range  of  variables  when 
using  Eq.  (8).  In  addition,  Eq.  (3)  is  somewhat  unwieldy  to  use,  and  Eq.  (5) 
is  valid  only  for  AQ/At  >  1  (or  B  >  2). 


An  approximate  approach  would  have  been  to  asstime  that 


[l  -  (Pc/Pr)]‘/Z  -  (Pc/2Pr)] 


so  Eq.  (6)  integrates  to 


2B 

Ps  ■  (B  +  2) 


r  {b  +  2)p  l 

r^rr 


•  (B+2)i 


As  the  equivalent  orifice -to-nozzle  area  ratio  B-*-i,  the  orifice  area 
restricts  the  gas  flow  (like  a  throttling  orifice)  so  P£  never  reaches  Pp.  To 
avoid  this  condition,  B  is  usually  sized  >  2. 


In  Fig.  3,  nor.iimensional  solutions  of  Eq.  (8)  are  presented  for 
yQ  =  0.  86  to  1.0,  B  =  1.  0  to  10.  0,  and  t  s  0.  005  to  i.  0.  Good  correlation 
between  the  analysis  and  the  data  for  five  gases  {both  at  sea  level  and  in 
vacuum)  is  noted.  These  results  indicate  that  the  time  required  to  reach 
steady-state  operation  (P  /P  =  1)  is  significantly  reduced  by  designing  the 
jet  reaction  chamber  with  a  small  chamber  volume  (Vc)  and  a  high  (B  >  Z) 
orifice -to-nozzle  area  ratio. 

From  Eq.  (8),  the  rise  time,  or  the  time  required  for  the  pressure 

in  the  thrust  chamber  to  reach  steady  state  (P/P  ~  i),  is  expressed  by 

c  s 

\  «(-B/a1)fn[B4y0  +  B)J  (9) 

During  shutdown,  the  orifice  valve  is  closed  and  the  gas  accumulated 
in  the  chamber  discharges  through  the  nozzle.  The  instantaneous  chamber 
pressure  during  an  isentropic  shutdown  is  given  by 

dPc/de  =  Jy.2(p c.)(‘-'()«V(pc)(3Y-l)/2v  ,  (10) 

Integrating  Eq,  (10),  the  pressure  ratio  at  time  0  becomes 

pc/pci  =  [l  -  -  Y)©]2Y/a~V>  dD 

V.  '  *_ 

with  the  results  as  given  in  Fig.  4. 
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DIMENSIONLESS  TIME  CONSTANT 


If  the  decay  time  0^  is  defined  as  the  time  required  to  discharge  from 


'  ### 

P  to  P  ,  then  from  Eq.  (ii) 

8  » 


«d  =  !<Pa/Ps)<1'V,/2V-,3^-  ')  (12) 

It  can  be  seen  that  0^  is  a  function  of  the  gas  properties  and  of  Vc/A^ 
(chamber  characteristic  length). 

The  correlation  between  the  ideal  isentropic  decay  and  the  experi¬ 
mental  results  is  not  quite  as  good  for  the  "slow"  (low  sonic  speed)  gases, 
such  as  the  Freons,  as  for  the  "fast"  gases,  e.  g.  ,  or  NH^  These  fast 
gases  discharge  so  rapidly  that  the  time  in  which  heat  transfer  anc’  other 
non-isentropic  imperfect  gas  effects  can  occur  is  restricted.  The  behavior 
of  the  more  complex,  heavy  Freon  molecules  indicates  that  the  decay 
transient  for  these  gases  follows  a  polytropic  process. 

In  a  polytropic  expansion,  deviations  from  isentropic  perfect  gas 
behavior  are  accounted  for  by  adjusting  the  ratio  of  specific  heats  y  through 
a  polytropic  efficiency  q  to  obtain  a  polytropic  exponent  k.  This  relationship 
becween  k  and  y  in  terms  of  qp  is  given  by 

k  =  y/[y  +  qp(i  -  y)]  (13) 


In  a  vacuum,  P  /F  ussumeu  equal  to  iO~J. 


The  thermal  expansion  efficiency  t)j  is  the  square  root  of  actual  drop  in  gas 
temperature  to  the  ideal  temperature  change  and  is  related  to  polytropic 
efficiency  by 


Values  of  tj  were  determined  from  the  test  results,  and  the  corresponding 
values  of  k  and  rjj  were  calculated.  The  results  are  listed  in  Table  1  for  ;he 
gases  tested  (except  which  is  so  fast  that  its  decay  is  well  described  by  an 
isentropic  decay).  Better  correlation  is  obtained  by  use  of  the  polytropic 
exponents  in  the  decay  transient,  as  shown  in  Fig.  4. 

B.  SYSTEM  CHARACTERISTICS 

The  thrust  and  nozzle  flow  transients  are  readily  determined  from  the 
chamber  pressure  history.  The  dynamic  impulse  bit  size  (total  impulse  per 
pulse)  and  the  dynamic  gas  consumption  (gas  consumed  per  pulse)  are 
obtained  by  integrating  the  rise  transient,  the  steady-state  value  (if  >  8^), 
and  the  decay  transient  of  thrust  and  of  gas  flow.  The  thrust  correlation^ 
factors  and  nozzle  discharge  coefficients  listed  in  Table  1  (and  discussed 
under  Performance  Losses)  were  used  in  the  transient  analysis  to  calculate 
impulse  bit  size  and  gas  consumption,  as  functions  of  command  pulse  ®c 
(defined  as  the  time  the  valve  is  open)  for  the  fire  gases  shown  in  Fig*.  5  r 
and  6.  Referring  to  Fig.  5,  the  difference  between  the  theoretical  and  the 
measured  varies  is  due  to  the  assumption  of  an  isentropic  decay.  It  was 
found  (Ref.  3)  that  the  correlation  is  improved  if  a  polytropic  expansion 


*  40  p*io  0t » 0.0188  in.  MEASURED _ 

*0.183  in*  Do ‘0.040  in.  oNHj  af'REON-14 


efficiency  of  about  85%  is  assumed  when  calculating  impulse  bit  size  and 
amount  of  gas  consumed.  However,  specific  impulse  is  not  affected  by  the 
assumption  of  isentropic  decay,  and  the  polytropic  approach  is  unnecessary 
for  this  calculation.  Note  that  in  the  region  of  small  command  pulses 
(8c  <  8r),  both  impulse  and  consumption  are  nonlinear  functions  of  command 
pulse.  This  nonlinearity,  or  deviation  from  the  ideal  square  pulse -wave 
(instantaneous  rise  and  decay),  is  produced  by  the  rise  and  decay  transients 
Equations  (8)  and  (11)  provide  the  designer  with  relationships  between 
the  variables  so  that  design  tradeoffs  may  be  evaluated.  For  minimum  gas 
consumption  at  any  given  command  pulse,  the  transients  should  be  mini¬ 
mized  or  possibly  designed  to  effectively  complement  each  other;  i.  e,  ,  the 
ideal  square  pulse  wave  is  approached  if  the  integral  of  the  chamber  pressure 
over  the  decay  time  can  be  designed  to  compensate  effectively  for  the  initial 
difference  between  the  product  of  the  steady-state  chamber  pressure  and  the 
rise  time  and  the  integral  of  the  chamber  pressure  over  the  rise  time,  i.  e.  , 


P  dO  should  approach  P  0 


-P 

r  vo 


P  dO  for  an  ideal  square  pulse  wave. 


The  dynamic  (pulsed)  specific  impulse  for  five  gases  as  a  function  of 
®c  (Fig.  7)  was  obtained  by  dividing  the  impulse  bit  by  the  weight  of  gas  con¬ 
sumed.  For  vacuum  operation,  specific  impulse  is  essentially  independent 
of  command  pulse  deration  ^  (unlike  I^and  w)  and  is  close  to  the  steady-state 
value.  Figure  7  indicates  good  correlation  between  the  analytical  and  experi- 


. 4.1 

xnciucu  udw. 

^For  command  pulses  less  than  20  msec,  the  data  becomes  nonrepeatabie 
because  of  valve  dynamic  limitations  and  increased  experimental  error. 
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C.  PERFORMANCE  LOSSES 


The  tests  in  this  study  were  made  using  a  15-degree  half-angle  conical 

aluminum  nozzle  with  a  measured  throat  diameter  Dt  of  0.  0188  in.  and 

an  expansion  ratio  of  48  to  1.  The  steady-state  chamber  pressure  was 

regulated  to  40  psia;  the  gas  temperature  (downstream  of  the  regulator)  was 

70°F*  The  ambient  pressure  in  the  vacuum  chamber  ranged  from  300  to 

1000  p,  depending  on  the  thruster  duty  cycle. 

The  perfect  gas  isentropic  flow  relationships  wej  e  used  to  calculate 

the  ideal  performance  (e.  g.  ,  characteristic  velocity,  nozzle  exit  pressure, 

and  temperature)  listed  in  Table  1.  The  tabulated  ideal  specific  impulse 

does  not  include  the  effects  of  nozzle  divergence  and  of  the  vacuum 

chamber  pressure.  Also  listed  in  1  sble  1  are  the  nozzle  throat  Reynolds 

number  (Nft  from  10,  000  to  53,  000),  the  nozzle  exit  Reynolds  number 

(Nre  from  3,  000  to  18,  000),  and  the  nozzle  exit  Knudsen  number  (N^  from 
-3  -3 

0.4  X  10  to  3.  3  X  10”  ).  Note  that  all  of  the  gases  operate  in  the  coninuum 

-E 

flow  regime  (N^  <  10  )  and  that  boundary  layer  viscous  effects  are  not  par- 

4 

ticularly  severe  (Nrt  >  10  ).  These  effects  may  become  important,  however, 
for  operation  with  values  of  N^  and  N  considerably  beyond  the  limits 
indicated. 

The  nozzle  performance  may  be  characterized  by  the  discharge 
coefficient  (the  ratio  of  actual  weight  flow  to  ideal  isentropic  weight  flow) 
and  the  exhaust  velocity  coefficient  (the  ratio  of  the  average  effective 
exhaust  velocity  to  the  ideal  isentropic  exhaust  velocity).  The  ratio  of  the 
actual  thrust  coefficient  to  the  ideal  thrust  coefficient  defines  the  thrust 
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correlation  factor  C^,  or 

C„  =  CJCt  ideal  =  C  C  . 
t  f  f  v  d 

It  follows  that  Cv  is  equal  to  the  impulse  efficiency,  or 


C  =  C.  /C  .  =  I/I.  ,  . 

v  t  d  ideal 


(16) 


The  values  of  C,  and  C,  in  Table  1  are  based  on  thrust  and  flow 
t  d 

measurements  taken  in  the  Aerospace  Guidance  and  Control  Laboratory. 

The  discharge  coefficient  of  the  test  nozzle  averaged  about  0.  85  for  all  the 
gases  except  NH^  and  Freon- 12.  The  relatively  high  of  nearly  1.  00  for 
these  two  gases  was  attributed  to  the  presence  of  fluid  condensate  in  the 
nozzle  flow.  The  thrust  correlation  factor  was  also  slightly  higher  for  the 
condensable  NH^  and  Freon- 12  gases.  The  average  impulse  efficiency  (Cv) 
was  about  91%.  The  efficiency  of  Freon- 12  is  slightly  lower  than  average. 
This  is  attributed  to:  heavy  molecular  weight,  uncertainties  in  thermo¬ 
dynamic  properties,  and  very  long  transients.  Values  of  C v  for  NH^  ranged 
from  88  to  95%,  depending  on  the  amount  of  condensation. 

A  method  of  correlating  the  losses  associated  with  small  nozzles  in 
terms  of  throat  Reynolds  number  and  nozzle  expansion  ratio  (Ref.  6)  was 


tt 


This  is  also  the  average  quoted  in  Refs. 


4  and  5. 
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investigated.  No  correlation  could  be  obtained  using  the  present  data  (a 
similar  lack  of  correlation  of  C£  and  C  in  terms  of  Reynolds  number  was 
reported  in  Ref.  7).  The  results  indicate  that  the  difference  between  tile 
ideal  isentropic  thrust  coefficient  and  the  measured  thrust  coefficient  is  a 
constant  of  about  0.  39  {except  for  NH^)  and  is  independent  of  the  throat  . 
Reynolds  number. 

In  addition,  an  analysis  of  the  boundary  layer  (Ref.  8)  was  made  to 
determine  its  effect  on  nozzle  performance.  The  effective  nozzle  expansion 
ratios  were  found  to  be  up  to  20%  less  than  the  geometric  ratios  (see  Table  1) 
with  a  consequent  performance  loss  of  about  3%.  It  was  estimated  that 
viscous  losses  and  changes  in  the  free -stream  flow  conditions  introduced  by 
the  boundary  layer  would  result  in  about  3%  additional  losses.-  The  total 
losses,  including  divergence  and  vacuum  backpressure  (but  not  condensation), 
are  estimated  to  be  about  9%. 

'  T  '  '  - 

*;  . .  I 

Isentropic  reaction  jet  expansions  of  NK^  and  Freon  -  .12  are  shbwn 
on  Mollier  diagrams  in  Figs,  8  and  9.  The  ideal  expansion  end-pointof' 
either  gas  is  well  into  the  two-phase  (liquid -vapor)  mixture  region  so  that 
a  metastable  supersaturated  condition  would  be  expected  to  occur  in  the 
nozzle.  The  NHj  expansion  not  only  crosses  the  saturation  line  but  also 
crosses  the  triple  line  (solid-liquid- vapor)  into  a  solid-vapor  region.. 

Since  incipient  condensation  of  steam  is  known  to  occur  at  the  5% 

5 The  data  in  Ref.  6  was  obtained  for  only  one  gas.  Hg,  and  at  Reynolds 
numbers  generally  less  than  10,  000,  whereas  the  values  of  Nrfc  for  the  five 
gases  tested  in  the  present  investigation  ranged' from  10,  000  to  53,  000. 
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Fig.  9.  Reaction  Jet  Cycle  with  Freon-  12 
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liquid  or  Wilson  line  (Ref.  9),  it  is  reasonable  to  assume  that  highly 
super-saturated  conditions  in  the  test  gases  cannot  be  sustained  without 
condensation  occurring,  thus  resulting  in  significant  deviation  from  an  ideal 
constant-entropy  expansion  process. 

The  problem  of  small  nozzle  condensation  was  first  observed  in  the 
present  study  when  NH^  and  Freon-12  were  tested  at  sea  level  (Ref.  1). 
The  existence  of  liquid  condensate  in  the  nozzle  discharge  was  demonstrated 
experimentally  (Ref.  3)  through  hot-wire  conductivity  measurements  and  by 


photographing  the  fluii  impingement  pattern  produced  on  a  black  matte  target 
downstream  of  the  nozzle  (Fig.  iO).  Attempts  at  direct  observation  of  the 
condensate  droplets  in  a  transparent  nozzle  through  the  Tyndall  effect 
(backscattering  of  a  light  beam)  were  inconclusive,  probably  because  the 
droplets  in  the  nozzle  were  smaller  than  the  wavelength  of  the  light  and 
were  moving  at  a  high  velocity  (s?rort  residence  time).  The  small  size  of 
the  nozzle  also  hampered  direct  optical  viewing  of  the  condensate. 

The  present  theoretical  treatment  of  nucleation  and  droplet  growth 
during  condensation  from  vapor  (i.  e. ,  the  disappearance  of  vapor  molecules 


and  the  formation  and  growth  of  fluid  droplets)  appears  unable  to  produce 
definitive  quantitative  results.  One  major  difficulty  concerns  nucleation 
kinetics  cr  the  reactions  (and  rates)  by  which  vapor  molecules  combine  to 
form  droplets  (Refs.  10  -  16).  The  treatment  of  the  influence  of  foreign 
impurities,  e.  g. ,  dust  or  ions,  which  provide  heterogeneous  nuclei  and 
promote  condensation,  presents  additional  difficulties.  In  the  present  study, 
the  pulsed  or  dynamic  mode  of  operation  introduces  another  transient  con¬ 
sideration.  Even  for  steady-state  operation,  the  small  nozzle  size  and  the 


short  residence time  (Ref.  15}  involved  (see  -Table .  i)  present  theoretical 

pe.  *■  *  ■  ■  ■  ■  .  - 

prohl*ih«ccmcerning  scale  effects.  l.f' -  ' 

^though  theproblem  of  determining  the  effect  Of  flow  condensation  on 
nozzle  performance  has  been  recognised  for  40  year 6  {Ref.  17),  the  experi¬ 
mental  work  inhomogeneous  nucleation  and  condensation  is  apparently  meager, 
particularly  with  respect  to  pulsed  flow  in  very  small  nozzles  *  Most,  of  the 
previous'  experimental  work  has  been  done  using .  saturated  water  in  large 
nozzles  (Refs.  9»  and  lb  -  20)  or  in  cloud  chambers  (Ref.  21).  The  qualita¬ 
tive  indications  are  that  the  . actual  flow  is  somewhat  greater  than  would  be 
theoretically  expected.  In  the  present  study,  the  observed  increase  in  Cd 
for  the  TIHj  and  Freon- 12  was  attributed  to  the  presence  of  fluid  condensate. 

Theoretical  methods  (Refs.  22  and  23)  for  calculating  the  condensation 
losses  were  examined.  These  results  indicate  from  3  to  7%  additional  per¬ 
formance  loss  for  two-phase  expansion  due  to:  energy  transfer  between  vapor 
and  droplets  (a  part  of  which  is  lost  by  dispersion);  the  appearance  of  super¬ 
cooled  states  and  impact  (shock)  condensation;  and- (possibly)  lack  of  time 
for  interphase  thermal  equilibrium. 

In  summary,  because  the  effects  of  condensing  vapors  on  small  nozzle 
expansion  processes  are  presently  not  well  understood,  additional 

"♦Accurate  measurements  are  usually  difficult  to  obtain  with  small  nozzles 
(Df  <  100  mils),  e.  g. ,  measurement  of  may  be  in  error  by  10%.  Also, 
microscopic  differences  in  construction  could  alter  the  nozzle  flow 
characteristics. 
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study  is  required  in  order  to  establish  accurate  means  of  performance  analysis 
and  characterization. 

D.  TRANSIENT  PERFORMANCE  APPROXIMATIONS 

When  the  rise  and  decay  thrust  and  flow  rate  transients  are  integrated, 
it  is  possible  to  define  time  integration  factors  so  that 
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s  d 


,•  -  (18) 


Using  these  empirical  factors,  dynamic  impulse  bit  size,  dynamic  gas  com- 

*  -  ‘  r  t  ^ 

sumption,  and  pulsea  specific  impulse  can  be  obtained  from  the  following 
approximate  relationships. 


Impulse  3it  Size 


For  0c  S  6r, 


W  *  CfWVr  +  KdV<W 


For  P  >  6  , 
c  r 


V=W.K“rt(,c-V,KA] 


(19) 


Gas  u.onsumption 


For  8  5  6  » 
c  r 


w  =  k  ajp  (k  e  +  K.e ,)  {9  /e  ) 

nts  rr  d  d  c  r' 


For  0  »  0  , 
c  r 


w  =  k  a  p  Tk  e  +  (0  -  e  )  +  k,9 ,1 

n  *  H  r  r  c  d  dj 


where 


K  =  choked  nozzle  flow  factor 
n 

--  a*  [2/(Y  + 


For  the  gases  tested  in  this  study,  values  of  -  0.  69  and 
=  0,  0685  were  found  to  provide  good  correlation  with  the  measured  data. 
Comparison  of  these  approximate  results  with  those  obtained  with  the 
isentropic  analysis,  the  polytropic  analysis,  and  the  data  are  shown  in 
Figs.  5-7. 

Specific  Impulse 

The  specific  impulse  is  the  ratio  of  Eq.  (19)  to  Eq,  (20),  or  simply 


/c^de 

TTaTI^  “  K 
*'  n  t  c  n 


=  0.  13  w&3  found  to  provide  good  correlation  with  an  isentropic  decay, 
wherers  K,  =  0.  0683  correlated  with  the  polytropic  decay  process. 


It  is  noted  from  Eq.  (21)  that  specific  impulse  is  theoretically  constant  and 


independent  of  command  pulse  0c> 


E. 


PERFORMANCE  COMPARISON  . 


Although  specific  impulse  is  the  most  important  single  measure  of  pro¬ 
pulsive  performance,  other  factors  become  important  for  system  design* 

As  pointed  out  in  Ref.  24,  consideration  should  be  given  to  the  propellant 
density  and  associated  tankage  weight.  This  may  be  accomplished  by  defin¬ 
ing  effective  system  specific  impulse  as  the  ratio  of  delivered  total  impulse 
to  the  combined  propellant  and  minimum  spherical  tankage  weight.  Thus 


2eff  =  I/[i  +  1.  5vP(p/<r)tJ 


(22) 


„  <\C  '  v. 


where 

X  - 

P  =  the  propellant  storage  pressure  (typically  3500  psia)~ 
v  =  the  propellant  specific  volume  v  'j 

-v'  '  o-, 

-  5  --v 

(p/«r)t  =  the  ratio  of  tank  density  to  working  stress  =  3.  25  X  10  ft  for 
titanium  with  a  safety  factor  of  2.  -  -fe? 

This  technique  was  applied  to  the  measured  .values  of  specific  UMpc&se 
for  the  five  gases  tested.  The  results  are  listed  in  Table  1.  The  condensable 
gasett,  NH^  and  Freon- 12,  were  a««umed  stored  as  liquids  at  their  room 
temperature  vapor  pressures  of  130  and  85  psia,  respectively,  and  Eq,  (18) 
was  modified  vo  account  for  minimum  gage  tanks,  based  on  actual  spacecraft 
design.  It- is  noted  that  eu  the  basis  of  specific  impulse,  has  the  highest 
performance  (260  sec),  whereas  on  the  basis  of  effective  system  specific 
impul&s,'  H^-has  She  lowest  pe^^-rmanceJ13  sec).  On  the  basis  of  effective 

the  highest  performance  (75  sec)  and  Freon-12 
is  at:2§  Mty&tpv&x.i  a  weight  penalty  for  the  heat  of 


r«nw.arMo‘ 


vaporization  of  NH^  and  Freon- 12  was  not  included  in  the  Table  1  perform¬ 
ance  values,  (In  the  laboratory  this  heat  was  provided  by  the  heat  capacity 
of  the  apparatus. )  In  the  case  of  NHy  about  500  Btu/lb  are  required  for 
vaporization.  This  heat  would  either  have  to  be  supplied  from  the  heat 
capacity  of  the  spacecraft  or  from  an  external  source,  depending  on  the 
application.  The  final  selection  of  one  propellant  over  another  would  depend 
on  the  results  of  a  detailed  design  study  for  a  specific  mission. 
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HI.  EXPERIMENTAL  PROGRAM  /• 


The  experimental  apparatus  shown  in  Fig.  1  includes.  air^gi^atei|’cold^ 
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gas  pneumatic  system,  a  thrust  stand,  and  an  electronic  command  ahd-._  „;-,p 

computer  system.  '  »  r.  ’;.X: 

... 

*  *  A,  .'ii1 

The  -old-gas  pneumatic  system  consists  of.  a  pressure 

,  < 

pellant  storage,  a  0  to  400  psi  pressure  gage,  a  thermocouple  tehap^.r4tttr#7if 

;  !•'  :■■■  . 

monitor,  a  pressure  regulator  adjustable  from  0  to  10Q  psi,  and  a. 
sure  monitor  gage.  When  the  system 


n  is  operajedbn  VLt‘  N-,  or 

a  •'  ‘  rv.  '  - 


pressure  and  temperature  measurements  <j*re  pied,  to  determine 

•V  5,  5 


ma  s  s  consumption  frog* 
system  is  operated  with  anhydrous  NIL*  or  F reon-tZ^the:  storage 

accurately  weighed  before  and  after  a  test  to  determine  the.maps 

. ’*  ■ 


vs  V 


'-S&tSksr 


during  the  tes*:. .  A  typical  test  consists  of  -  a e  vegal  thousand  '  jgre 
command  pulses.  •  " 

r  *-  v  .  „  -  f--VA  ;-:<a  .  -o* o^bisss^ssSk 

*  J  . '  ■  ’  Is  .JR;  . 

_  *  fSssHi. 


The  thrust  measurement  apparatus.  Fig.  .1  lv^has-A;  noi« 

*'•  Jsk  « .  :  s-tatognt 

the  end  of  a  short  length  of  tubing  to  serve  the  dual  fun^tioiirdf ‘i.‘pneui^ttcX;j;:i. 

.  .  ‘  ‘'“-.."a;  V  V  -  •»? 

feed  line  and  a  cantilever  beam  spring*?;  When  thrust  in.prodUcp,ffi,v.rftie  fcjSsime.: 
is  displaced  in  an  amount,  proportional  to  the  thrust  mag^4bi<4^^3^  d 
of  the  beam  is  sensed  fcy  a  position  transducer,,  and  an 


Xi 

portional  to  the  position  (thrust  level)  is  generated, 
is  provided  through  the  shearing  action  of  a  0a|  plate 

fluid.  A  high,  natural  frequency  of  the  cantilever  beam  and  nozzle  assembly 

<  /-  C  ' 

(250  cps)  is  necessary  to  provide  adequate  dynamic  response  for  the  thrust 
transients. 
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Fig.  1L  Thrust  Stand  Assembly 


The  natural  frequency  of  the  thrust  measurement  fixture,  was  verified 
experimentally  by  two  methods.  The  first  was  displacement  of  the  beam  by 

\  •  -  i 

a  discharge  of  gas  from  the  pneumatic  system  through  the  nozzle.  The  natural 
frequencies  of  the  fixture,  both  with  and  without  damping  fluid,  v^ere  displayed 
on  the  oscilloscope  and  photographed.  A  photograph  of  the  undamped  oscilla¬ 


tion  is  shown 


noid  valve  in  sue 


»n  in  Fife  42»  -  The.  8ec.qnd.me.thpdi^aa  ppjiUqnipg  a  coaxial  *ole- 


igjtatrge  through  the 


valve  and  impingefib^lhie''tl^^  fi^tdi«  hoxaf*  blocfei.  ;  Tf*£:£e  suiting  displace - 

<Vt  -  *'  -r.  '#* 

ment  of  the  fixturi^ajr  ^fefiSMpi&j&oiK^^  adjusting  the  high 

_ -  ;  —  _ SMI, V-*-  -  *'■  S-i  :-i - .a. - - 


pressure  air  suppw(^i^tyeg-  was  set 

equal  to  the  , 

3B£L*&  Aii- '  K  ’•  5?&  • 

13  is  a  photograph  o4  the  lraha»e»t  ros|^^t^t%  ihrawt^j^nd  excited  Hjp- 

siV^  .  r’t  V.  ♦-'•»'  -V.i r-V  ».  V  <.?& 

such  an  external  force.  The  thrust  stand  dynamic  response  is  more  thin 


adequate  for  the  rise  times  encountered,  as  evidenced  by  the  cor  relation  of  . 

,  ‘  v  £***'■’  J'"~  '< 

thrust  and  pressure  traces  in  Fig.  14.  (A  pressure  transducer  with  a  >JK 

l  .  - 

1200  cps  acoustic  frequency  was  used  to  monitor  pressure, )  „  c  ■  <-* 

.  4a>t*  i 

The  thrust  stand  was  calibrated  in  two' ways:  (1)  using  a  Sch *jc$S&kH\' 

,  -  *  {*V''  sgj  \ 

%  *  «J  Vi  ' 

Tumico  force  gauge  and  providing  a  scale  factor  in  terms  of  volts, pe r  pound 

of  force;  (2)  orienting  the  thrust  stand  so  that  the  beam  was  subjected  to  a 

'  'M  ...  ,  - 

1  g  loading.  Test  weights  were  then  placed  on  the  cantilever  beam.anda 
correlation  of  force  versus  electrical- output  (deflection)  was  made.  Both 

-  !r  A  v 

methods  were  in  excellent  agreement.  v 

‘  •  c 

A  schematic  of  the  command  and  computer  system  is  shown  in  Fig.  15. 


A  preset  counter  or  clock  was  used  to  fix  the  repetition  rate  of  the  command 


ELECTRICAL 
COMMAND  PULSE 
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Fig.  14.  Correlation  of  Thrust  and  Chamber  Pressure 
for  Hydrogen,  (steady-state  thrust  =  0.0146  lb 
steady-state  chamber  pressure  =  40  psia; 
command  pulse  =  50  msec) 
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pulses,  and  a  digital  counter  recorded  the  number  of  cycles  or  pulses  applied 
to  the  solenoid  valve.  The  output  of  the  thrust  transducer  was  amplified, 
demodulated,  and  fed  to  the  voltage -to -frequency  converter.  The  output  of 
the  converter  was  fed  into  a  reversible  counter  which  was  used  to  integrate 
the  thrust  profile  ana  thus  obtain  the  total  impulse  per  pulse  (impulse  bit  size). 
The  reversible  counter  computed  the  impulse  bit  including  drift  or  offset  dur¬ 
ing  a  time  increment  At.  Then,  the  counter  totalled  the  integrated  offset  and 
drift  for  an  equal  time  increment  At  and  obtained  the  difference  between  the 
two  values  to  provide  a  direct  readout  of  total  impulse.  The  reversible 

'k  •  V  ’ 
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counter  provided  improved  accuracy  over  the  two  separate  counters  pre¬ 
viously  used  (Ref.  3)  for  the  same  function  by  eliminating  any  errors  due  to 
differences  in  timing  sensitivity  between  individual  counters.  '  . 

The  specific  impulse  was  obtained  by  dividing  the  total  impulse 'by  the 
gas  consumed  as  measured  over  a  large  number  of  command,  pulses  ,(>1000). 
The  solenoid  valve  closing  time  was  adjusted  to  equal  the  valve  opening  time 
by  valve  driver  circuitry  modifications  (Ref.  29).  Therefore,  this  configura¬ 
tion  resulted  in  a  valve  open  time  equal  to  the  electrical  command  pulse  width, 
with  a  3.  5-msec  delay  relative  to  the  command  pulse.  ,  - 

To  complete  the  experimental  evaluation,  the  repeatability  of  impulse 
bit  size  was  verified.  Two  thousand  20-msec  command  pulses  were  applied 
to  the  solenoid  valve,  and  the  total  impulse  of  each  pulse  was  recorded.  The 
Hi BtrtHntfnn  nin<-  of  th*»s**  data.  1A.  indicate?  the  excellent  repeatability 

of  impulse  bit  size.  -  • 

An  estimated  error  analysis  of  the  experimental  iconfiguration  was  made. 
This  analysis  included  the  effects  of  mass  measurement  errors,  limited 
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dynamic  response  of  the  thrust  fixture,  calibration  errors,  and  electronic 
errors.  The  error  estimate  took  into  account  the  variation  in  experimental 
measurement  accuracy  with  command  pulse  width  and  with  the  propellant  • 
used.  The  average  error  for  all  propellants  and  command  pulse  widths  is 
approximately  5%,  the  typical  spread  of  measured  data  points.  Table  1 
presents  the  estimated  measurement  error  associated  with  each  propellant 
for  command  pulse  durations  of  20  and  200  msec.  The  predicted  measurement 
accuracy  of  5%  correlates  well  with  the  experimental  data.  In  Fig.  ?  for 
example,  80%  of  the  experimentally  measured  data  lies  within  5%  of  the 
theoretical  predictions. 

During  the  course  of  the  experiments,  a  leakage  problem  became 
apparent  when  testing  H^.  Because  of  the  low  molecular  weight  of  H-,  (see 
Table  1),  a  leak  of  only  1.  25  X  10“^  lb  of  gas  while  obtaining  a  data  point 
would  produce  a  20%  error  in  measured  specific  impulse.  The  same  quantity 
of  leaking  N2  would  result  in  n  5%  error  in  specific  impulse.  It  was  found 
that  the  leak  was  caused  by  the  O-ring  used  to  seal  the  nozzle  in  the  nozzle 
block  and  could  have  been  prevented  by  careful  installation. 

An  additional  problem  was  encountered  during  the  H2  experiments. 
Significant  variations  in  specific  impulse  were  observed  during  the  initial 
tests.  Subsequent  investigation  revealed  that  these  variations  were  due  to 
the  effects  of  gaseous  contamination  occurring  taring  the  filling  of  the  pro¬ 
pellant  tank.  Normally,  when  the  propellants  were  changed  (».  g.  >  changing 
from  N2  to  H,,),  the  complete  pneumatic  system  was  evacuated  using  a" 
vacuum  pump  to  remove  residual  gases...  Then  the  propellant  tank  was  filled 


with  the  new  gas.  However,  ambient  air  in  the  4-ft  long  fill  hose  had  rot 
been  considered.  Calculations  indicated  that  the  error  in  s  ecific  impulse 
introduced  by  the  contamination  of  with  trapped  air  could  be  as  large  as 
20%  due  to  the  change  in  the  average  gas  molecular  weight.  This  conclusion 
was  experimenta’Iy  /e.ified.  The  propellant  loading  procedure  was  then 
modified  to  include  a  vacuum  purge  of  the  propellant  fill  line  as  we  1  as  of 
the  pneumatic  system. 


IV.  CONCLUSIONS 

Accurate  predictions  of  impulse  bit  size,  gas  consumption,  and  effective 
specific  impulse  can  be  made  using  the  gas  c»/namical  relationships  given  for 
the  transient  pressure  histories.  These  relationships  constitute  an  effective 
analytical  tool  for  optimizing  system  performance  l^r  a  given  control 
nirement. 

Both  the  gas  consumed  and  the  impulse  bit  size  are  nonlinear  functions 
,  command  pulses  for  small  pulses  but  become  linear  when  steady-state 
c.hamr,  uiessure  is  reached.  The  nonlinearity,  or  deviation  from  the  ideal 
square  pulse-wave,  is  due  to  the  rise  and  decay  transients,  For  minimum 
gas  consumption  at  a  given  command  pulse,  the  transients  can  be  minimized 
by  proper  design  so  that  the  rise  transient  effects  are  essentially  offset  by 
the  decay  transient  effects.  A  simplified  approach  to  transient  performance 
determination  given  by  Eqs.  (19)  through  (21)  can  be  used  for  preliminary 
design  purposes. 

The  pulsed  vacuum  specific  impulse  is  essentially  independent  of 
command  pulse  width  and  is  close  to  the  steady-state  value,  A  loss  of  f'-om 
5  to  10%  in  specific  impulse  was  observed  for  command  pulses  less  than 
20  msec.  This  degradation  is  attributed  to  the  increasing  dominance  of 
solenoid  valve  dynamics.  Because  of  inci easing  interest  in  small  pulse- 
width  performance  (<20  msec),  additional  study  in  this  area  is  required  in 
order  to  characterize  solenoid  valve  dynamic  effects  on  performance. 


Ofl 


This  discharge  coefficient  of  the  nozzle  tested  averaged  about  85%  for 
all  the  gases  except  NH^  and  Freon-1?-.  The  average  impulse  efficiency  was 
about  91%.  No  correlation  of  the  nozzle  losses  could  be  obtained  in  terms  of 
Reynolds  number.  The  results  indicate  a  constant  loss  in  ideal  isentropic 
thrust  coefficient  of  approximately  0.  39  (except  for  condensing  gases,  MH^ 
and  Freon-12).  Additional  study  is  warranted  on  the  subject  of  small  nozzle 
expansion  processes  and  on  the  performance  effects  of  condensing  vapors. 

On  the  basis  of  specific  impulse  alone,  shows  the  highest  perform¬ 
ance.  If  propellant  and  tank  weights  are  considered,  then  NH^  has  the  highest 
effective  system  specific  impulse.  Since  the  NH^  was  assumed  to  be  stored 
as  a  liquid,  the  effect  of  the  required  heat  of  vaporization  would  have  to  be 
considered  in  spacecraft  design  applications.  The  final  selection  of  pro¬ 
pellant  should  be  based  on  a  detailed  study  of  the  specific  mission  require¬ 


ments. 
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13  ABSTRACT 

The  pulsed  propuisi/e  performance  of  low-thrust  reaction  jets,  typical  of 
those  used  for  small  spacecraft  attitude  control,  is  analyzed  and  compared 
with  the  results  of  laboratory  experiments.  Five  gases,  hydrogen,  nitrogen- 
ammonia,  Freon-12,  and  Freon-14,  are  investigated  using  a  48  to  l  expansion 
ratio  nozzle.  The  transient  processes  which  dominate  the  short-pulse'  or 
limit-cycle  mode  of  thruster  operation  are  formulated.  These  relationships 
show  good  correlation  with  the  data.  The  apparatus,  procedures,  and 
techniques  required  to  obtain  accurate  test  results  for  a  low-thrust,  dynamic 
mode  of  operation  are  described.  Impulse  bit  size,  gas  consumption,  and  ~ 
specific  impulse  are  characterized  in  terms  o £  thruster  geometry,  gas 
properties,  and  command  pulse  width  to  provide  a  basis  for  dptimum  system 
design.  A  simplified  method  for  calculating  dynamic  impulse  bit  size, 
dynamic  gas  consumption,  and  pulsed  specific  impulse  as  a  function  of 
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the  gases  tested  is  evaluated  u/  a  technique  which  includes  the  influence  of 
tank  and  propellant  weights,  as  well  as  specific  impulse. 
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